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I. INTRODUCTION 



The production of particles in the interacting nuclear medium has been a rich field of 
study The particle production is a quantity which can be to a large extent estimated 
in theoretical models. The particle production was used as a tool of the study of the highly 
nonequilibrium dynamics in the early stage of the nuclear collision and also as a probe of the 
behavior and a characteristic of the almost equilibrated system reached in the course of the 
collision. The estimation of the particle production probability usually is based on the free 
production cross section with some possible medium modifications included. This production 
cross section, taken for the 2-nucleon collision process, is used in the semiclassical Boltzmann- 
Uehling-Uhlenbeck (BUU) type simulations of the nuclear dynamics. A completely different 
approach to the particle production emerges from the quantum transport equations. The 
importance of the quantum transport description was realized in the correct treatment of 
the Landau- Pomeranchuk-Migdal effect in medium . The approach based on the quantum 
transport equations takes into account the finite time extent between successive collisions 
and regularizes the divergent cross section for long-wavelength photons which require large 
formation time. 

Two problems related to the particle production in intermediate energy heavy ion col- 
lisions are particularly interesting : the subthreshold meson production and the medium 
modifications of the meson properties. The production of the particles below the nucleon- 
nucleon energy threshold requires the presence of some additional mechanism in nuclear 
medium allowing for the production. Besides the obvious Fermi motion of nucleons, ad- 
ditional mechanisms could include the initial high momentum component in the nucleus 
and the many nucleon collisions. These phenomena can be treated in a unified way in the 
quantum nonequilibrium Green's function formalism 0-^. The medium modification of the 
meson properties include the possible mass shift and width modifications in medium. If this 
occurs, a natural question arises also about possible strong medium modifications of the 
production rates. In the case of very strong medium modification of the meson properties 
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the meson production should be treated in a quantum transport framework. 

In this work it will appear that the difference between the quantum and semiclassical 
production rates is very important. Usually one estimates the off-shell effects by taking 
a generalized parameterization of the fermion spectral functions in the calculation of the 
meson production rate in 2-nucleon collision. However, in the quantum transport equation 
the perturbative expansion using resumed fermion propagators is organized differently. The 
first nonzero contribution to the meson production comes from the one-loop diagram and 
its contribution is generally large. Another way of discussion of the subthreshold particle 
production is to estimate the particle production in 3 or more nucleon collisions processes. 
A large class of such meson production processes involving many nucleons is included in the 
one-loop meson self-energy with interacting off-shell fermion propagators. 

In the present work in order to make a detailed estimation of the quantum and semi- 
classical production rates in the same system we will use a simple model. In this model 
of self- interacting nonrelativistic fermions weakly coupled to mesons we will calculate the 
one-loop quantum meson production rate and the corresponding semiclassical rate from the 
2-nucleon collision process. Clearly we cannot use the free on shell meson production cross 
section as a starting point. Instead we will use a simple meson-nucleon coupling and we 
will derive the corresponding semiclassical production rate. The numerical calculations of 
the dynamics of the fermion Green's functions are done in a similar way as in the work 
of Danielewicz 0. The production rate of mesons is estimated in a perturbative way, i.e. 
we assume that it is small and that the meson production does not modify the fermion 
dynamics. 

In order to make the calculation feasible we assume that the system is spatially homoge- 
neous. The quantum evolution equations are much simpler in this case, since the momentum 
integrals are the same as in the usual evaluation of perturbative diagrams and only the time 
integration must be done on a contour in the complex time plane. From the experimental 
evidence it seems as well established that the particle production around the threshold, but 
sufficiently above the absolute threshold, is an incoherent process. The number of produced 
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mesons scales very well with the number of participants. The result of the calculation in the 
homogeneous system giving the particle production rate per unit volume can be scaled to 
the size of the actual system under consideration. Thus, we believe that our results could 
be used for actual estimates of the production probability. Of course this would require the 
use of realistic meson nucleon vertices, and also the inclusion of the delta degrees of free- 
dom. Our goal in this work is restricted to the presentation of the particle production in the 
nonequilibrium transport equations, with comparison to the semiclassical production rates. 
However, as a good motivation one should keep in mind the fact that the comparison of the 
particle production estimated in the off-shell Green's functions framework would represent a 
powerful test of the quantum effects. The particle production is one of the rare observables 
which can be reasonably well estimated in a homogeneous system and then compared to 
the experimental production rate per participant. Only the estimation of the absorption 
of the produced mesons must be done in a finite system. The difference in the transport 
coefficients between the quantum and the semiclassical evolutions discussed up to now |^-|^ 
is not directly observable and could be effectively included in a medium modification of the 
nucleon-nucleon cross section. 



II. PARTICLE PRODUCTION IN QUANTUM TRANSPORT EQUATIONS 

We will consider a system of nonrelativistic interacting fermions. The first numerical 
solution of the Kadanoff-Baym equations for a homogeneous system was presented in Ref. 
. For the dynamics of the fermions we follow the same approach. In sections |IIB| and [II C 



we will calculate the meson production rate in the system of interacting off-shell fermions. 



A. Quantum transport equations for fermions 

In this section we recall the notation of |^,^ and the numerical procedure for the solution 
of the nonequilibrium dynamics of fermion Green's functions. The Hamiltonian of fermions 
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interacting with a two-body potential V, can be written using the annihilation and creation 
operators ijj^x) and as follows : 

H = J dx ^^(x) ^ ^(x) + - J dx J dy ^\x)tjj\y)V{x — y)xl){y)^{x) . (1) 

The spin and isospin indices are not written explicitly and we take a scalar potential for 
simplicity. The one-body evolution of the system can be described by the Green's functions : 

G^{xi,ti,X2,t2) = i < i'Hix2,t2)^Hixi,ti) > , 

^^(Xi, ti, X2, t2) = -i< ^HiXi,ti)^lj{x2, t2) > , (2) 

where the field operators tpH are in the Heisenberg representation and the average < . . . > 
is the average over the initial state. These Green's functions are a particular case of the 
general contour time ordered Green's function : 

G{xi, ti,X2, t2) = i < Tc1pHix2, t2)lpHixi,ti) > , (3) 

where the field operators are ordered according to the time on the contour in the complex 
plain (Fig|T]). For ti, t2 on different branches of the real time part of the contour the Green's 
function G can be replaced by two Green's functions G*^*^^^ defined on the real line. In a 
spatial homogeneous system there is no dependence on the macroscopic coordinate ^^^^^ , 
and after Fourier transforming in the relative coordinate xi — X2 we obtain the Green's 
function in momentum space ti, ^2) • The evolution of the system starting at some 

time to can be described by the Schwinger-Dyson equations : 

d [ill 
iTr-G{p,ti,t2) - ujpG{p,ti,t2) = 6c{ti,t2) + / dtJ:{p,ti,t)G{p,t ,12) 
oti Jc 

-t^Gip, tiM) - ^pGip, ti, ta) = Scih, ts) + / dtG{p, t )^{p, t M) , (4) 
012 Jc 

where ^1,^2) is the self-energy which must be calculated in some approximation, 
'^0(^15^2) is the delta function on the contour and the time integration t' is done on the 
contour (Fig. [l|). We neglect the mean field potential and the medium modification of the 
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nucleon effective mass, so that = For ^i? ^2 real and on different branches of the 
contour we obtain the Kadanoff-Baym equations : 

Oil \Jto+iTo Jto / 

+ r diT.<{p,tiJ)G<{p,t M) 



+ / + / )dtj:<{p,ti,t)G>{p,t,t2) 

\Jto Jt2 J 

i-^G>{pMM) - ujpG>{pMM) = ( r + r)dtT>{pM,t )G^{p,tM) 

oil \Jto+iTo Jto J 



+ / dt'S>(p,ti,t')G>(p,t',t2) 

/ /-to -JTO ftQ\ 

^ L )^^^^(P'^i'0G'>(p,t,t2) , (5) 

and similar equations for the evolution in ^2- The description of the fermion evolution using 
only the one-body fermion Green's function involves the assumption that there are no 2 or 
many nucleon correlations present at the initial time U0[3 . Such many- nucleon correlations 



could be important for some processes involving nucleons in nuclei [|TT . 

The actual solution requires first the solution on the imaginary time part of the contour, 
which is done iteratively. Then the Kadanoff-Baym equations are solved on the real time 
part of the contour. Also similar equations are solved for the case where one of the time 
arguments of the Green's function is on the real part and the other one is on the imaginary 
part of the contour. Of course the Hamiltonian of the evolution on the imaginary part of 
the contour must be modified. Otherwise we would always obtain the ground state of the 
one-body evolution with the Hamiltonian (|l]), for tq 00. This can be done as an additional 
constraint or some modification of the Hamiltonian on the imaginary part of the contour. 
We will describe the nuclear collision as the equilibration of two counter-streaming flows of 
nuclear matter. Accordingly we choose, following 0, as the initial state two Fermi spheres 
of radius pj = 255MeV/c centered at ±EcniL ^jth the kinetic energy of the Hamiltonian 
on the imaginary part of the contour Pt~^(\pl\^j>cou/2) jj^g^g^^^ £_ ^^le real part of the 
contour. The resulting evolution of the fermion Green's function is the so called evolution 
with correlated initial state. The self energy is taken in the direct Born approximation : 



G<>\p2, h, h)G<>\p + pi - P2, t2, h)V{p2 - pf , (6) 

which allows for an efficient numerical evaluation. The factor 4 comes from spin isospin 
trace and the potential V{p) = vr^/^?7^Voexp(— ^) (Vq = AbSMeV , rj = 0.57 fm) should be 
thought of as a local approximation to the T matrix. 

Another possibility is to use only the real part of contour for the calculation of the 
evolution of the system. The fermion momentum distribution at t = to is taken as two 
Fermi spheres. The kinetic energy in this case is increased in the course of the evolution and 
the final distribution corresponds the equilibrium distribution with higher temperature than 
in the case of the evolution including also the imaginary part of the contour. This real time 
only quantum transport evolution of the fermion system will be refered to as the evolution 
from an uncorrelated or a Hartree initial state. 

We have used an axially symmetric system for the evolution of the Kadanoff-Baym 
equations, with momentum grid of 46MeV/c and of the size \p\ < 920MeV/c. The imaginary 
time evolution was done for tq = 3/m/c and the real time evolution was done for 20/m/c. 
The momentum integrals are done by a Fourier-Bessel transformation and the time integrals 
are done in a usual numerical integration for each pair of times ti, t2 (ti > ^2)- 

B. Meson production 

The meson production can be analyzed in a model where in addition to the Hamiltonian 
H for interacting fermions (|l|) we have a meson-fermion interaction vertex and of course the 
meson kinetic energy term. For the meson-fermion interaction we could take the form of the 
interaction as in the pion-nucleon case [|1^ : 

Hi = g Jdx(^tp(x fa''di(l)j{x)T^tp(x)+^p(x)^a^di4>j{x)W^^p(^^ , (7) 

where is the pion field operator, r and a are the isospin and spin matrices. The interaction 
strength g incorporates in general a momentum dependent form factor. However, in our 
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study the detailed structure of the interaction vertex is not important. Both in the one- loop 
diagram and in the semiclassical production rate the meson self-energy acquires the same 
coefficient from the meson-fermion interaction vertex. In the case of the interaction vertex 



in the momentum space, where q is the momentum of the meson. The interaction denoted in 
the following by could be used also for any kind of meson-nucleon vertex, after taking the 
spin isospin trace. In the present study, which is a first estimation of the quantum production 
rate, we will not be interested in the detailed predictions. Instead we will calculate the 
rates of the meson production scaled by 1/ for the comparison of the quantum transport 
equations and the semiclassical equations. Also we will treat the meson mass as a free 
parameter, performing the calculation for several different masses of produced particles. 

The evolution of the nucleon meson system can be analyzed in the same way as the 
system of fermions, discussed in the previous section. The meson Green's functions can be 
defined as : 



where in the following we restrict our study to one meson species. The meson nucleon inter- 
action induces a backreaction of the produced mesons on the fermion evolution. However we 
will suppose that the meson-nucleon vertex is weak, and that the meson production does not 
modify the fermion dynamics. The fermions do interact by the interaction V which induces 
the equilibration of fermions and leads to nonzero width of the fermion Green's functions. 
Thus, the meson production can be calculated from the known fermion Green's functions. 
We can write also the Kadanoff-Baym evolution equations for mesons : 



written above (|^, we have a common factor : 



(8) 





(9) 



1 




ti, h) - %D{p, ti, t2) = Sc{ti,t2) + I dtD{p, ti, t')n(p, t,t2) , (10) 
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where 11 is the meson self-energy and Vtp = ^ ml + is the relativistic meson kinetic energy. 
We take at this point the meson kinetic energy in the relativistic form in order to have a 
correct estimate of the energy of the meson for any meson mass and momentum. Otherwise 
the mesons are treated in an analogous way as the nonrelativistic fermions. The meson self- 



energy is taken in the one- loop form (Fig. |^). The equations (10) can be used for a complete 
study of the meson evolution in the nonequilibrium nuclear matter, including the in medium 
modifications of the meson properties, the meson production and absorption. In the present 
work we are interested only in the calculation of the meson production rate. It is equal to 
the rate of the change of the equal time meson Green's function iD'^{p, t, t) = dN{p, t)/d^p, 
which is the meson momentum distribution in a homogeneous system . Supposing that the 
number of produced mesons is small, which is true for sufficiently small A^, we can neglect 
the meson occupancy on the right hand side of the equations (PH]), taking : 

D<{p,h,t2) = 

Do>(p,ti,t2) = -^exp(-^^]p(^l -ta)) . (11) 

This approximation is consistent with the neglect of the meson influence on the fermion 
dynamics. The production rate for mesons takes the form 



dN{p, t) I rt 

d^pdt 



27^e|^-^^ dtU<{p,t,t)D^{p,t ,t) 

pto — iro I I , \ 

+ rftn<(p,t,t)Do>(p,t,t) . (12) 



The meson one-loop self energy is : 

n<(p,ti,t2) = -^A2/ ^G<(p-g,ti,t2)G>(g,t2,ti) • (13) 

At this point, some comment is again welcomed concerning the integration on the imag- 
inary part of the contour in the Kadanoff-Baym equations. We assume that the imaginary 
part of the evolution does not create mesons. This means again that we modify the Hamil- 
tonian on the imaginary part of the evolution on the time contour, or equivalently that the 
evolution operator on the imaginary part of the contour includes a projection on the states 

9 



without mesons |TJ]. In this way we have no mesons at the initial time to. The problem of the 
presence of mesons in the initial state is by its own very interesting and could be addressed 
by the Kadanoff-Baym equations (^) and (p!0|). The calculation of the meson distribution 
in the initial state would require a self-consistent solution of the interacting meson-fermion 
system. This is beyond the scope of the present work, where we restrict ourselves to the 
perturbative calculation of the meson production rate from the fermion Green's functions. 

The rate of the meson production (|12D consists of two time integrals, a real time and an 
imaginary time integration. For t ^ oo the imaginary time integration is negligible, and 
we obtain as the production rate the imaginary part of the meson self-energy taken at the 
energy Qg. For finite integration times, the situation is different. In particular, taking the 
stationary fermion Green's functions on-shell : 

G'<(p,ti,t2) = if{p)exp{-iujp(ti -t2)) , 

G>{pMM) = f{p))QW{-i^p{ti-t2)) , (14) 

we obtain from the real part of the time integration in Eq. (|T2p an oscillatory behavior : 

dN{g,t) + (15) 



d?qdt J {2nY — ujg-p + Up 

which only on average gives a zero production rate. On the other hand, taking the same 
form of the on-shell fermion Green's functions ( |14|) also on the imaginary interval of the 
time integration in Eq. (|12D and taking tq ^ oo one obtains : 

(-) 

exactly for each t. Thus, we believe that the imaginary part of the time integration is 
important in the calculation of the production rate in a finite time interval after to- Of 
course in the system that we solve numerically we have off-shell Green's function which 
gives a nonzero production rate. Also in the far from equilibrium system at the beginning 
of the collision the Green's function depend strongly on the macroscopic time . One 
should remember that even in the stationary case the analytical continuation of the self- 
energy from the imaginary to the real times is nontrivial and in fact it is possible only in 
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the perturbative calculation, when the analytical expressions for the self energy are known. 
The method presented above is way of calculating the real time quantities from the Green's 
functions evolved in the imaginary time. It would be interesting to analyze in details this 
problem in the case of a finite temperature stationary solution, i.e. in the case where the 
imaginary evolution is the usual finite temperature evolution over an imaginary time interval 
1/T with the same Hamiltonian on the real and the imaginary part of the contour. 

C. Numerical results for the quantum production rate 

The numerical calculation of the fermion evolution was performed for the energy of the 
collision of ISOMeV and llOMeV/nucleon in the cm.. We have calculated two evolutions 
of the fermion Green's functions, corresponding to a correlated or uncorrelated initial state. 
The evolution was followed up to 20/m/c. At that time the fermion distribution is almost 
spherical and equilibrated. Also, at t = 20/m/c there is already no influence of the initial 
part of the time contour on the evolution. However, the evolution for the initial correlated 
and the evolution for the initial Hartree state have different temperatures at large times. 
The temperature in the second case is larger and so is the meson production at large times. 
The calculation of the one-loop meson self energy is done in a similar way as for the fermion 
self-energy, i.e. by a Fourier-Bessel transform calculation of the convolution momentum 
integral in Eq. ([T3|) . For zero meson momentum the integration is done directly. 

We have calculated the meson production in three ways. The first possibility is to take 
the Hartree initial state and the whole evolution of the fermion Green's functions is in the real 
times. For this case we calculate the meson production using Eq. ([12|) with only real time 
integration on the right hand side. For the case of the correlated initial state the evolution 
of the fermion Green's functions is performed both on the real and imaginary part of the 
contour in Fig. ^ The corresponding meson production rate can be calculated in two ways. 
One way is to integrate only over real times on the right hand side of Eq. (|1^), the other one 
is to take the whole expression (^). Of course, at large times, the two approaches give the 
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same result, since the contribution of the imaginary time integration is neghgible for large 
times. However, at the beginning of the evolution, the differences are important. In most 



of the cases, the full expression (|1^) is well behaving, with less oscillatory behavior. The 
oscillatory behavior at the beginning of the evolution is due to the finite time integration. 
This gives only on average an expression for the production rate. 

In Figs. 1^, H and |^ is shown the time dependence of the scaled production rate p-^^i^ 
for three different masses of the mesons 20, 140 and SOOMeV respectively. We prefer to 
use the scaled production rate, since it allows us a general discussion. In particular, we can 
analyze the behavior around g = where the production rate would be if multiplied by the 
pion vertex strength. The general behavior is as expected, starting with oscillations until 
the asymptotic production rate settles in. This initial time dependence of the production 
rate is not due to the nonequilibrium character of the evolving system. Rather, it is due 
to the finite time integration in the expression for the production rate, as discussed above. 
For not very large values of the meson mass or meson momentum q, i.e. not very far 
from the mass shell, we can parameterize the fermion Green's functions using a Lorentzian 
spectral function. In that case one expects that the oscillations in the production rate would 
be damped with a time l/F, where F is the average width of the spectral function. In our 
case the width is around 50MeV which should give a decay time of Afm/c in accordance 
with the Figs. |], ^ and |^. The asymptotic value of the production rate at large time is not 
zero. However, for large mass of the meson or large momentum of the meson the asymptotic 
value of the production rate is much smaller then the amplitude of the oscillations at the 
initial time. In order to get stable results for the case of large mass and/or momentum 
of the produced meson we had to use very small time step {.067 fm/c and .15/m/c for 
the real and imaginary evolution respectively). This allows to get a correct estimate of 
the fermion spectral function far from the mass-shell. The asymptotic production rate of 
mesons is systematically larger for the evolution from the uncorrelated initial state (dashed 
lines in Figs. ^, ^ and |^). It is a manifestation of the the higher temperature reached in the 
equilibration of the uncorrelated initial state. In Fig. |^ is shown the time dependence of the 
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meson production rate for large times. For the case of meson of mass 140 and 300MeV an 
almost stationary production rate is reached. However for the mesons of mass 500MeV still 
oscillations of high amplitude are seen. These oscillations are slowly damped to a positive 
value corresponding to the asymptotic meson production rate. 

In Fig. 1^ is shown the production rate of mesons at large time (t = 20/m/c), as a 
function of the transverse momentum. At this time the production rate of mesons is almost 
isotropic. The quantum meson production rate for large meson energy was estimated as the 
average over its last oscillation in the time dependent meson production rate. Similarly as 
in the plots of the time development of the production rate, the production rate is larger 
for the uncorrelated initial state. One should note that the plots in Fig. |^ do not represent 
the genuine transverse momentum distribution of produced mesons, since the interaction 
strength can in general depend quite strongly on the momentum of the produced meson 
q. In Fig. ^ is shown the total number of mesons produced during the evolution as a function 
of the transverse momentum. The highest number of mesons is produced in the evolution 
from an uncorrelated initial state. In the time integrated production rate, a big difference 
appears between the meson production calculated from the full expression (|T2[) and the 
meson production calculated from the Eq. including only the real interval of the time 
integration. The difference can be observed in the plots of the time dependence of the meson 
production rate for small times. The large oscillations of the production rate at small times 
give after time integration a large contribution. This contribution is unphysical, since in 
a nuclear collision there is no sudden switching on of the interaction at time t = Iq. The 
time integrations are always infinite and the production rate should not have the oscillations 
characteristic of the finite time integration. At each time the locally asymptotic, but not 
necessarily stationary, meson production rate can be defined, if the time integration would 
be much larger than the inverse fermion single particle width l/F. In the nonequilibrium 
Green's evolution presented in this work, we have to start from a finite initial time. The 
procedure with imaginary time integration in the initial state leads to a better description 
of the evolution of fermions and of the meson production rate. 
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Close to the mass-shell the fermion Green's functions are often parameterized as : 

G<(p,ti,t2) = -^(l-/(p))exp(-zcup(ti-t2)-rp|ti-t2|) , (17) 

and one gets for the meson production rate : 
dN{q = Q,t) 2 f d^P 



,H,t ' lW^'^^'-^^'^''<fn- '''' 

This formula gives Lorentzian dependence on the meson mass if the width Tp is weakly 
depending on momentum. In Fig ^ is shown the dependence of the asymptotic production 
rate of mesons at zero momentum as a function of the meson mass. This dependence is 
close to an exponential and cannot be approximated by a Lorentzian. In the appendix |B] is 
shown that even if the fermion spectral function is Lorentzian we do not get a Lorentzian 
dependence of the meson production rate on the meson mass, if we use a more general 
expression for the equilibrium Green's functions than the Eq. (|1^. 

The fermion Green's functions obtained as a result of the numerical solution of the 
Kadanoff-Baym equations are very different from the parameterization (|T^ The function 
f{p) in the above parameterization corresponds to the momentum distribution of fermions. 
However, we have found that the imaginary part of -iG<(»(p,ti,t2) exp(ico'p(i(:i — ^2)) is 
nonzero and gives an important contribution to the meson production rate. The Lorentzian 
parameterization of the spectral function is expected to be valid near the mass shell. Ac- 
cordingly, it cannot be applied for the description of the meson production which involves 
the far from the mass-shell spectral function of fermions. 

III. COMPARISON TO THE SEMICLASSICAL PRODUCTION RATE 

The semiclassical meson production rate from a homogeneous BUU-like evolution of the 
collision will be presented. The BUU evolution is done for the Hartree initial momentum 
distribution and for the initial momentum distribution obtained from the imaginary time 
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evolution of the Kadanoff-Baym equations. Such calculations correspond to the standard 
approach to the meson production in nuclear collisions. 

Other results will also be presented for the semiclassical production rate with fermion 
momentum distribution taken from the equal time Green's functions generated in the quan- 
tum evolution. It is interesting to compare the semiclassical meson production rate in a 
quasistationary state to the quantum production rate in the same system. 

The semiclassical production rate is defined in App. 0. In this chapter we use the 
expression ([ATD and (|A2|) , with 7 = 50MeV. 



The semiclassical evolution equation for fermion densities can be obtained from the 
Kadanoff-Baym equations [@-|^. We take the nucleon-nucleon scattering cross section cor- 
responding to the direct Born fermion self energy : 



The initial momentum distribution of fermions f{p, to), should be taken as two Fermi spheres 
separated in momentum space by the momentum of the collision. However, we know that 
the nuclei in the ground state are not Hartree states with with zero temperature. In par- 
ticular they have a large momentum component in the fermion momentum distribution. In 
our model involving spatially homogeneous system, it can be described by taking as the 
initial state the Hartree state evolved on the imaginary time contour by the Kadanoff-Baym 
equations. This corresponds to allowing for the initial quantum correlations, which modify 
the momentum distribution. This effect is expected to be important for the subthreshold 
meson production rate. Accordingly, we will use also the fermion momentum distribution in 



A. Semiclassical evolution equations and meson production 




(27r)'^5^(pi +P2-P3- P4)27T6{uJp^ + UJp^ - Up^ - WpJ 



f{pi,t)f{p2,t){l-f{p,,t)){l-f{p,,t)) 




(19) 
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the correlated state as the initial momentum distribution to)- The average kinetic en- 
ergy in the correlated state is around lOMeV/nucleon larger than in the Hartree state. The 
momentum distribution in the correlated initial state has a large momentum component. 

The BUU-like evolution following the semiclassical fermion collision term (|19]) was solved 
on the same momentum grid as the quantum evolution. The time step was taken equal 
to 0.5/m/c. In the effectively 5 dimensional integration on the right hand of Eq. ([T9|) 3 
dimensions were integrated by Monte-Carlo integration and two momentum integration were 
summed over the momentum grid. 

From the semiclassical fermion momentum distribution obtained from Eq. (p^Of), we can 
calculate the semiclassical meson production rate using Eq. Another possibility is to 



use the fermion momentum distributions obtained from the quantum evolution equations : 

f{p,t) = -tG<ip,t,t) . (20) 

As the system in the quantum evolution reaches an almost stationary state after 20/m/c of 
real time evolution, we can compare the meson production rate in this stationary state as 
obtained from the quantum and the semiclassical expressions for the meson production. The 
two stationary states that we discuss below have different temperatures and correspond to 
the initial correlated state and the initial Hartree state in the quantum evolution equations. 



B. Discussion of the results for the semiclassical meson production rate 

The time dependence of the semiclassical production rate is very similar for different 
meson masses or momenta. Thus, we present only one Figure of the time dependence 
of the semiclassical production rate. In Fig. ^ is shown the time dependence of the 
semiclassical meson production rate for the meson momentum q = 184MeV/c and for three 
different meson masses. At the initial time the semiclassical production rate has two values, 
depending on the initial state of the evolution. For the correlated initial state (dotted and 
dashed-dotted lines in Fig. |T^) the production rate at the beginning of the evolution is 
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higher than for the Hartree initial state (sohd and dashed hnes in Fig. |T0|). It is due to 
the fact that the average kinetic energy in the correlated initial state is higher than in the 
Hartree state. 

The BUU-like evolution from the initial Hartree state, has the lowest average kinetic 
energy. For the higher energies of the produced mesons, especially for the case of Qg = 
352MeV (the lower panel), the semiclassical production rate increases during the evolution. 
Such a behavior is expected, since the system approaches thermal equilibrium with some 
small thermal high momentum tail. However, as the average kinetic energy is conserved, 
the final equilibrium state has the lowest semiclassical meson production rate among all the 
results presented in Fig. 0. 

On the other hand, the quantum transport evolution starting from the initial Hartree 
state, leads to an increase of the kinetic energy. The initial correlation energy is changed 
into the kinetic energy As a result the average kinetic energy is increased by 

SOMel^/nucleon and it is the largest among all the different evolutions considered in Fig. 
|IU|. Thus, it is obvious that also the semiclassical production rate would be the largest in 
this case. Within the first few fm/c the semiclassical production rate increases very fast. 
After this initial increase the meson production rate increases only slowly. This slow change 
of the semiclassical meson production rate is determined by the equilibration of the momen- 
tum distribution of fermions in the later stage of the Kadanoff-Baym equations evolution of 
fermions. 

The two evolutions starting from the correlated initial state have very similar average 
kinetic energies during the whole evolution. Thus, the semiclassical meson production rate 
for the BUU-like evolution and for the quantum evolution from a correlated initial state are 



not very different (the dotted and dashed-dotted lines in Fig. [T^). The BUU-like evolution 
starting from the correlated initial state can be considered as a model of the BUU evolution 
of the initial nucleon momentum distributions in nuclei that includes the large momentum 
tail of quantum origin (dotted line in Fig. |r^). This kind of calculation represents all 
the quantum effects that can be taken into account in the approach using the 2-nucleon 
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process only and on-shell particles. To go beyond the discussion of the quantum effects 
in the initial state only, i.e. to include the off-shellness of the fermion Green's functions, 
requires an off-shell calculation of the meson production rate as was presented in Sect. p^lB 



and 111 Q The difference between the BUU evolution from the Hartree initial state and the 



BUU evolution from the correlated initial state is especially important for the case of the 
production of meson of mass SOOMeV. This meson mass corresponds to the subthreshold 
meson production, since the energy of the produced meson is 352MeV and the average cm. 
energy is ISOMeV (the kinematical limit for the meson production in the Hartree initial 
state is around 450MeV). The ratio of the semiclassical meson production in the correlated 
initial state evolution and in the Hartree initial state evolution changes from 26 at t = to 



3 at t = 20/m/c . From Fig. |T0] we see that the BUU-like evolution gives similar results 
as the the quantum evolution if the correlated initial state is used. The use of the initial 
correlated state is justified in the Kadanoff-Baym equations. It is a much more realistic way 
of starting the quantum transport equations evolution than the Hartree initial state. 

We do not show in Fig. ITUI the quantum meson production rate since the quantum rate 



is one order of magnitude larger. Moreover, the asymptotic quantum meson production rate 
can be defined only for large times. The large time quantum and semiclassical meson pro- 
duction rate as a function of the transverse momentum of the produced meson are compared 
in Fig. 0. The quantum meson production rate is about one order of magnitude larger 
than the semiclassical meson production rate. The dashed and the dotted line represent the 
quantum and semiclassical meson production rate for the larger temperature state reached 
in the evolution from the Hartree initial state. The solid and the dashed-dotted line repre- 



sent the same for the evolution from the correlated initial state. From Fig. 11 we see that 



the semiclassical and the quantum meson production rates are very different for the same 



state reached in the evolution of the system. This is very different from Fig. |T0|, where 
we have seen that the differences between the semiclassical and quantum evolutions of the 
correlated initial state are not big if we are using only the semiclassical meson production 
rate formula. 



In Fig. |T2| we present the transverse momentum distribution of the meson production 
rate at t = 20/m/c for the lower energy of the coUision (llOMel^/nucleon in the cm.). The 
three "physically" well founded approaches are the quantum evolution from the correlated 
initial state, the BUU evolution from the correlated initial state [|T2| and the BUU evolution 
from the Hartree initial state (the solid, dashed-dotted and dotted lines in Fig. O). Besides 
these results we plot the semiclassical meson production rate for the fermion momentum 
distribution calculated in the quantum evolution of the correlated initial state. Similarly 
as for the higher energy the semiclassical meson production rate is similar for the BUU 
and quantum evolved correlated initial state. Also similarly as in the results presented in 



Fig. |TT| the quantum meson production rate is about one order of magnitude larger than 
the semiclassical meson production rate. The difference is even larger if compared to the 
semiclassical meson production rate for the BUU evolution of the Hartree state. We have 
used the same effective width 7 = 50MeV in the calculation of the matrix element for 
semiclassical meson production, as for the higher energy, since the single particle fermion 
width is reduced only by a few percent, when changing the initial energy. We do not show the 
results for = SOOMeV, because no asymptotic meson production rate could be estimated 
for this meson mass and large meson momentum. 



IV. CONCLUSIONS 

We have presented the first approach to the calculation of the meson production rate in 
the non-equilibrium quantum transport equations. The Kadanoff-Baym equations for the 
evolution of fermions were solved and the resulting off-shell fermion Green's functions were 
used for the calculation of the one-loop meson self-energy. The one-loop meson production 
rate shows strong oscillations at the beginning of the evolution of the system. This behavior 
comes from the fact that the time integral in the meson evolution equation is finite. Only 
after some time the meson production rate reaches some asymptotic value, which can be 
defined as the quantum meson production rate. We have found that the behavior of the 
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meson production rate at small times improves if the imaginary part of the time integration 
in the meson evolution equations is taken into account. We plan to study this question in 
details for the case where the fermion Green's functions on the imaginary and real parts 
of the time contour is the finite temperature Green's function corresponding to the same 
stationary finite temperature state. 

The formalism based on the off-shell nucleon propagators allows to estimate the effects 
of the initial quantum correlations on the meson production rate. We have found that these 
initial correlations are important and that their inclusion increases the semiclassical meson 
production rate. This effect is especially important for the subthreshold meson production 



{iris = SOOMeV in Fig. |TOD. The initial correlations in the fermion distribution can be 
evolved with BUU-like fermion evolution equation. Also the fermion momentum distribu- 
tion can be obtained from the equal time fermion Green's functions calculated with the 
Kadanoff-Baym equations. Both calculations give similar average fermion kinetic energy 
and similar semiclassical meson production rate. Thus, it makes no dramatic difference in 
the semiclassical particle production if the nucleon momentum distribution is evolved ac- 
cording to the Kadanoff-Baym equations or according to the BUU equations, if the initial 
momentum distribution includes the quantum correlations. 

The situation is very different if the quantum meson production rate is used to calcu- 
late the number of produced particles. In the system reached by the equilibration of the 
nuclear collision of the energy ISOMeV^/nucleon in cm, the quantum meson production 
rate is usually one order of magnitude larger then the semiclassical meson production rate. 
This is an important difference which requires detailed attention when discussing the par- 
ticle production in a hot nuclear matter. The phenomenological consequences for the pion 
production may be less important since we discuss only the direct pion production pro- 
cess. The direct pion production is known to be less important than the delta production, 
at least in the semiclassical approach. The delta production should be less influenced by 
the nucleon off-shellness. However, as a result of our estimates the direct pion production 
should not be neglected in the total pion production yield if the quantum production rate 

20 



is used for the direct pion production. More detailed calculations using a model incorpo- 
rating delta, nucleon and pion degrees of freedom is required to address quantitatively the 
phenomenological predictions for pion production. More work is also needed in the study of 
the parameterization of the retarded nucleon propagator which comes into the calculation 
of the semiclassical meson production rate. This could be studied first in the equilibrium 
state using the quantum transport equations. 

In view of our results, the models of the particle production in the hot nuclear medium 
should be based on the Kadanoff-Baym equations, or should be based on off-shell transport 



equations ||l6l. Great care must be shown in the use of approximate fermion spectral func- 
tions. The numerical calculations show that the simple parameterizations of the fermion 
spectral function as in Eq. [B^ does not give the correct quantitative result. This is even 
more pronounced for the less general parameterization of the fermion Green's function in 
Eq. ([17D , giving a qualitatively different dependence on the meson mass. 
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APPENDIX A: SEMICLASSICAL PRODUCTION RATES 

The semiclassical production rate can be obtained from the cut diagrams, by putting 
the Greens functions G^^'^^ on shell flj5|jr5[|. This procedure leads to zero for the one-loop 
diagram (Fig. 0) that was discussed for the quantum production rate. The 3-body energy 
conserving delta function has no solution for our kinetic conditions. In the next order in 
V, only a vertex correction to the one-loop diagram appears. This also does not give any 
contribution on-shell. In the second order in the interaction V there appear 8 direct dia- 
grams (Figs. O, andlT^). There are 8 other exchange diagrams, that we will not discuss 



here since for the fermion evolution only direct diagrams are taken. The diagrams in Fig. 

21 



T3| should not appear in the perturbative expansion, since they represent self-energy correc- 
tions on fermion lines. These are already included in the dynamically calculated fermion 
Green's functions. However, they should appear naturally in a perturbative expansion of 
the production rate in vacuum, using the on-shell Green's functions. The separation of the 
3-body event which is the scattering with meson production from the succession of a 2-body 
scattering and successive meson emission can be done by the non-zero width in the retarded 
propagator in the diagrams in Figs. 0, |l^ and |l^ [l^. Of course the 2-body scattering 



and successive meson emission gives zero contribution on shell. If such a procedure is used, 



the diagrams in Fig. [T^ are a first term in a series of diagrams involving the direct Born 
self-energy insertions on fermion lines. The sum of the diagrams in such a series should cor- 
respond to the one-loop meson-self energy diagram with resumed fermion Green's functions. 
One diagram of the next order in this series is shown in Fig. 16. The cutting procedure on 
the higher diagrams in this series leads to the processes involving many fermions on-shell. 
Such many nucleon processes should be taken into account in the calculation of the sub- 
threshold particle production. The many nucleon processes, which can be obtained from the 
cutting of the direct Born self energy insertions, are resumed in the off-shell Green's func- 
tions calculated in Sect, p.1 A|. In the following we will take only the diagrams in Fig. Ofor 



the calculation of the semiclassical meson production rate. This correspond to taking only 
the 2-nucleon term in the expansion of the meson self energy in the number on interacting 
nucleons. The corresponding matrix elements is given below (j^^) (Fig. |18D. The diagrams 



in Fig. |l^ are zero in the case of the interaction vertex and also more generally for a 
large class of theories. 

Identifying each cut propagator with the fermion Green's function on shell {iG^^^^ de- 
pending on the direction of crossing the cut), and the propagator on the left (right) side of 
the cut with the retarded (advanced) propagator we obtain the following form of the meson 
production rate : 

^5^ = 7 (2^7 (2^7 (2^7 (2^"^l (2vr)5(p,+p,-p3-P4-g) 
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27r5{up, + Up, - ujp, - - n,)f{p,,t)f{p2,t) (1 - f{ps,t)) (1 - /(p4,t)) . (Al) 

Here analogously as for the calculation of the quantum production rate (sect. [11 B[ ) we put 
the meson Green's function on shell and neglect the meson occupancy in the calculation of 
the production rate. The square of the matrix element |Mp in the above expression can be 
identified with the sum of the cut meson self-energy diagrams (Fig. |TB|, |T1| and |T^) with the 
cut fermion Green's functions iG'^'^^^ factored out. The square of the matrix element |Mp 



can be written in the form of the diagram in Fig. [T^ where the fermion Green's function 



is the retarded propagator and the external Green's functions have momenta corresponding 



to the integration momenta in Eq. ( [Al|) . The diagrams in Fig. |Tj and |T5] appear also in 
the expansion of the self-energy in the quantum transport equations, i.e. in the expansion 
using the resumed fermion propagators, but at a higher order in V . Thus in order to com- 
pare the semiclassical rate to the one-loop quantum rate, only the cut self-energy insertion 
diagrams (Fig. 0) should be taken. Clearly, even using the finite width for the retarded 
propagator in the semiclassical calculation no correspondence is found between the order of 
the perturbative expansion on shell and off shell. The cut diagrams in Fig. Ill and ITB| give 



contributions of the same order as the cut self-energy insertion diagram (Fig. |T3|). However, 
in the quantum expansion theses diagrams are of different order. Moreover, already at the 
first order in a vertex correction to the one-loop quantum diagram appears, which is still 
zero on- shell. 

For the diagrams in Fig. |1^ we obtain (Fig |T8|) : 

1 



^ Up^^g + Up2 - Up.^ - Up, + 27 

+ 1 ^ -?] ■ (A2) 

Up^ + + Up^+q - Up, +1-1 ) 



After the calculation of the traces over spin and isospin indices, the interaction g'^ and a 
factor appear. This can be written in the same way as for the one-loop calculation, 
i.e. as A^. The additional factor 4 comes from the spin isospin trace in the second fermion 
loop. In most of the cases, unless otherwise stated, we use the production rate with constant 
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width parameter 7 = 50MeV. This is the typical average fermion width as obtained from the 
quantum transport equations for fermions at the colhsion energies considered here. Of course 
in general it is momentum and time dependent, but in the calculation of the semiclassical 
production cross section we take a constant value of the fermion width. 



APPENDIX B: MESON PRODUCTION IN EQUILIBRIUM 

In this section we will present a simple estimate of the meson production in an equilibrium 
fermion system at finite temperature. For the fermion spectral function we will take the 
Lorentzian form. The fermion Green's functions, after Fourier transforming in the relative 
time variable, will be parameterized in the following way : 

where we take a constant width Tp and 

^»<-' - l + e.p((L,)/T) («^) 

is the equilibrium Fermi distribution at finite temperature (T = 30MeV, n = AOMeV, 
r = 60MeV, \p\ < gOOMeV^, |w - tOp] < ISOOMeV). 

We have calculated the meson production from the one loop diagram : 

The above equation is different from the Eq. ([T^) used in Sect. P C| . Here we are using the 



general equilibrium form of the Green's functions 0] with a Lorentzian spectral function. 
As a result, the dependence of the meson production rate on the meson mass is almost 
exponential (see Fig. |T^, similarly as in the numerical results in Fig. P. This dependence 



is not Lorentzian as expected from Eq. ([T8|). Even for Lorentzian spectral function with 



constant width the resulting dependence of the production rate on the meson mass is not 
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Lorentzian. This makes impossible a direct estimation of the single particle fermion width 
from the comparison of the meson production rate for different meson masses. 



In Fig. |T9| is also shown the semiclassical meson production rate calculated from the Eq. 

l1|) at zero meson momentum. For the equilibrium fermion momentum distribution we 

have taken the equilibrium momentum distribution at finite temperature : 

^\ 
2m 



/(p) = /oh?-) (B4) 



or the momentum distribution including the quantum effect 

/(p) = -^/^G<(p,a;). (B5) 

In the second case the rate of the meson production is much larger, reflecting the large 
momentum tail in the momentum distribution calculated from Eq. (p5|). Similar large 
momentum tail is present in nuclear matter at zero temperature. This effect was found to be 
important for the semiclassical meson production in our nuclear collision model (Sect. |11 C| ). 
In Fig. |T^ are shown also the semiclassical meson production rates for the width parameter 
7 — i> 0. This strongly modifies the meson production rate for small meson mass. The meson 
production rate has a singularity for rris 0. However, the nonzero 7 parameter changes 
this behavior. The meson production rate stays finite for rris if 7 > 0. This is similar to 
the Landau-Pomeranchuk-Migdal effect discussed in 0|. The semiclassical production cross 
section in medium is modified and has no singularity, even-though the meson production is 
given by the semiclassical expression 

The quantum meson production rate is much larger than the semiclassical production 
rate in Fig. p. However, such a comparison is not very meaningful. The Green's functions 
G^^^^ used in the one-loop meson self energy (|BT|), where not calculated selfconsistently 
for a system of fermions interacting with potential V. On the other hand the interaction 
potential V enters in the expression for the semiclassical production rate of mesons. 
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FIGURES 
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FIG. 1. The contour in the complex time plane on which the Green's functions are defined. 
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FIG. 2. The one-loop meson self-energy diagram. The solid lines are the fermion propagators, 
the dashed lines indicate the location of the meson-nucleon vertices. 
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FIG. 3. Quantum rate of the production of mesons of mass 20MeV as a function of time. The 
dashed hne is the meson production rate for the the uncorrelated initial state. The sohd hue is the 
meson production rate for the correlated initial state with both real and imaginary time integration 
in the meson production rate (12). The dashed-dotted line is the meson production rate for the 



correlated initial state with only real time integration in the meson production rate (12). The 
upper, middle and lower panel correspond to the transverse momentum of the produced meson of 
0, 184 and 414MeV/c respectively. The energy of the collision is ISOMeF/nucleon. 
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FIG. 4. Quantum rate of the production of mesons of mass lAOMeV as a function of time. 
Details are similar as in Fig. ^ 
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FIG. 5. Quantum rate of the production of mesons of mass SOOMeV as a function of time. 
Details are similar as in Fig. ^ 



31 




t [fm/c] 



FIG. 6. Quantum rate of the production of mesons at zero momentum as a function of time. 
The dashed hne is the meson production rate for the the uncorrelated initial state. The sohd hne 
is the meson production rate for the correlated initial state with both real and imaginary time 
integration in the meson production rate (12). The dashed-dotted line is the meson production 
rate for the correlated initial state with only real time integration in the meson production rate 
([l^). The upper, middle and lower panel correspond to the mass of the produced meson of 140, 
300 and bOOMeV respectively. The energy of the collision is ISOMeF/nucleon. 
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FIG. 7. Quantum rate of the production of mesons as a function of the transverse momentum 
at the end of the evolution (t = 20 fm/c). The dashed hne is the meson production rate for the 
uncorrected initial state. The solid line is the meson production rate for the correlated initial 
state. The upper, middle and lower panel correspond to the mass of the produced meson of 20, 
140 and SOOMeV respectively. The energy of the collision is ISOMeF/nucleon. 
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FIG. 8. The number of produced mesons as a function of the transverse momentum at the end 
of the evolution (t = 20/m/c) for the quantum one-loop meson production rate. The dashed line is 
the number of produced mesons for the the uncorrelated initial state. The solid line is the number 
of produced mesons for the correlated initial state with both real and imaginary time integration 
in the meson production rate (0). The dashed-dotted line is the number of produced mesons for 
the correlated initial state with only the real time integration in the meson production rate ([T^. 
The upper, middle and lower panel correspond to the mass of the produced meson of 20, 140 and 
300MeV respectively. The energy of the collision is ISOMey/nucleon. 
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FIG. 9. Quantum rate of the production of mesons as a function of the meson mass at the end 
of the evolution {t = 20/m/c), for mesons with zero momentum. The dashed hne is the meson 
production rate for the uncorrelated initial state. The solid line is the meson production rate for 
the correlated initial state. The energy of the collision is ISOMeF/nucleon. 
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FIG. 10. Semiclassical rate of the production of mesons at momentum q = 18AM eV/c as a 
function of time. The sohd and the dotted hnes represent the semiclassical meson production rate 
for the BUU-like evolution of the fermion distribution from a Hartree and a correlated initial state 
respectively. The dashed and the dashed-dotted lines represent the semiclassical meson production 
rate for the quantum transport equations evolution of the fermion distribution from a Hartree 
and a correlated initial state respectively. The upper, middle and lower panel correspond to the 
mass of the produced meson of 140, 300 and SOOMeV respectively. The energy of the collision is 
ISOMeF/nucleon. 
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FIG. 11. Rate of the production of mesons as a function of the transverse momentum at the 
end of the evolution {t = 20/m/c). The dashed hne is the quantum meson production rate for the 
uncorrelated initial state, the dotted line is the semiclassical meson production rate for the same 
final fermion momentum distribution. The solid line is the quantum meson production rate for the 
correlated initial state, the dashed-dotted line is the semiclassical meson production rate for the 
same final fermion momentum distribution. The upper, middle and lower panel correspond to the 
mass of the produced meson of 20, 140 and SOOMeV respectively. The energy of the collision is 
ISOMeF/nucleon. 
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FIG. 12. Rate of the production of mesons as a function of the transverse momentum at the 
end of the evolution {t = 20/m/c). The sohd line is the quantum meson production rate for the 
correlated initial state, the dashed-dotted line is the semiclassical meson production rate for the 
same final fermion momentum distribution. The dashed line is the semiclassical meson production 
rate for the BUU evolved correlated initial state. The dotted line is the semiclassical meson 
production rate for the BUU evolved Hartree initial state. The upper and lower panel correspond 
to the mass of the produced meson of 20 and lAOMeV respectively. The energy of the collision is 
llOMey/nucleon. 
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FIG. 13. The cut diagrams for the meson self-energy consisting of self-energy insertions on the 
fermion lines in the one-loop diagram. The wiggly line denotes the time local interaction, the cut 
is indicated by the long dashed line, other symbols as in Fig. 
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FIG. 14. The cut diagrams for the meson self-energy which are not self-energy insertions on 



the fermion lines in the one-loop diagram. Symbols as in Figs, y and 13. 
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FIG. 15. The cut diagrams for the meson self-energy which are not self-energy insertions on the 
fermion lines in the one-loop diagram. These diagrams are zero for the pion-nucleon interaction 



vertex. Symbols as in Figs, g and |13 



41 




FIG. 16. A cut diagram for the meson self-energy of a type of direct Born self energy insertions 
on fermion lines. Many other self energy insertion can be added, giving after cutting processes 
involving many on-shell fermions. Symbols as in Figs. § and 
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FIG. 18. The matrix element from the cut fermion self-energy insertion diagrams from Fig 
text for details. 
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FIG. 19. Rate of meson production at zero momentum as a function of the meson mass. The 
dashed Hne is the quantum one- loop result ( p3| ) . The dotted line is the semiclassical production rate 
with the quantum equilibrium fermion momentum distribution and the dashed-dotted line is the the 
semiclassical production rate with the semiclassical equilibrium fermion momentum distribution 
(7 — s- 0). The solid line close to the dotted line is the same as the dotted line but with a non-zero 
width (7 = 60MeV) of the retarded propagator in the calculation of the semiclassical production 
cross section. The solid line close to the dashed-dotted line is the same as the dashed-dotted line 
but with a non-zero width (7 = 60MeV) of the retarded propagator in the calculation of the 
semiclassical production cross section. The Fermi energy, the temperature and the width are 40, 
30 and 60MeV respectively. 
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